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1 INTRODUCTION 

The past fifteen years have witnessed some of the more memorable discoveries in the 
ninety-year-old field of superconductivity, including the high-Tc oxides in 1986 JIl, 
the alkali-doped fuUerenes in 1991 0, the charge-injected fuUerenes in 2000 0, and 
the binary compound MgB2 in January 2001 [Q. The discovery of superconductivity 
at such high temperatures (40 K) in the simple s,p-metal compound MgB2 was quite 
unexpected. The absence ||^ of the problematic weak-link behavior of the high-Tc 
oxides and the relative ease of synthesis in various forms p has raised hopes that 
MgB2 may be suitable for numerous technological applications. 

To aid in the search for related compounds with even better superconducting 
properties and to help identify the pairing mechanism, a great deal of current research 
is dedicated to fully characterizing MgB2 in both its normal and superconducting 
states. A wide range of experiments, including isotope effect R, M, heat capacity R 



rof , inelastic neutron scattering [0, 0, NMR |]T3|, and photoemission spectroscopy 
I^ , support the picture that MgB2 is a phonon- mediated BCS superconductor in 
the moderate coupling regime. The fact that the B isotope effect is fifteen times that 
for Mg H] is clear evidence that the superconducting pairing originates within the 



graphite-like B2-layers, consistent with electronic structure calculations [15, 16, 12 



TTj ?] whereby MgB2 is a quasi-2D material with strong covalent bonding within the 
boron layers. The anisotropy in the superconducting properties is appreciable, the 
upper critical field ratio H"2/H^2 reportedly being 1.7 [|l8l or 2.7 |19|, but far less than 



that observed in the high-Tc oxides A full characterization of all anisotropic 

properties awaits the synthesis of sufficiently large single crystals. 

High pressure studies traditionally play an important role in superconductivity. 
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Even without a detailed understanding of why changes with pressure, a large 
magnitude of the pressure derivative dTc/ dP is a good indication that higher values 
of Tc are possible at ambient pressure through chemical means. It is not widely 
appreciated, however, that the pressure dependence Tc(P), like the isotope effect, 
contains valuable information on the superconducting mechanism. In fact, in simple 
s,p-metal BCS superconductors like Al, In, Sn, or Pb, is found to invariably 
decrease with increasing hydrostatic pressure or isotopic mass; in both cases 
the reduction in Tc arises from changes in the lattice vibration spectrum, electronic 
properties having minimal effect. In transition-metal systems the isotope coefficient 
may deviate considerably from the BCS value a = 0.5 and the pressure dependence 
Tc{P) is determined by changes in both lattice vibration and electronic properties. 

Soon after the discovery of superconductivity in MgB2, three groups reported 
independently that Tc decreased under the application of high pressure, but the rate 
of decrease varied considerably. Lorenz et al. ||2^ carried out ac susceptibility 
measurements in a piston-cylinder cell to 1.8 GPa using the fluid pressure medium 
Fluorinert FC77 and obtained dT^/dP ~ -1.6 K/GPa. Saito et al. |]23[ measured the 
electrical resistivity p to 1.4 GPa using a similar pressure technique with Fluorinert 



FC70 and reported the pressure derivative -1.9 K/GPa. Monteverde et al. |24 
extended the pressure range to 25 GPa in resistivity measurements in an opposed 
anvil cell with solid steatite pressure medium; three of the four samples studied 
exhibited widely differing pressure dependences with initial values of dT^/dP ranging 
from -0.35 to -0.8 K/GPa. When pressure is applied to a solid pressure medium 
like steatite, the sample is subjected to sizeable shear stresses which may plastically 
deform a dense sample or compact a loosely sintered sample, as in the experiments of 
Monteverde et al. [^. Shear stresses are known to influence the pressure dependence 
of Tc, particularly in elastically anisotropic materials such as the high- Tc oxides |25 



or organic superconductors |]2B|. Fluid pressure media such as Fluorinert, methanol- 
ethanol or silicon oil remain fluid at RT to a certain pressure, but soon freeze upon 
cooling at temperatures well above the superconducting transition temperature Tc ~ 
40 K of MgB2, thus subjecting the sample to shear stresses, albeit relatively small 
ones. Only helium remains fluid at 40 K for pressures to 0.5 GPa. 

Since in electronic structure calculations Tc is determined by the unit cell dimen- 
sions and atom positions, for a quantitative comparison with theory it is essential to 
complement the determination of Tc(P) with accurate measurements of the pressure 
dependence of the structure parameters. Effects of pressure on the structure can 
affect the superconducting transition temperature through changes in the electronic 
structure, phonon frequencies, or electron-phonon coupling. For structures where the 
pressure effects are isotropic, changes in the electronic structure are usually subtle be- 
cause the Fermi energy and features of the Fermi surface tend to simply scale together 
with the cell volume. However, when the compression is anisotropic, as a result of 
significantly different bonding strengths in different crystallographic directions, large 
pressure-induced changes in the electronic structure can occur. For example, in the 
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layered copper oxide superconductor HgBa2Cu04+a;, where the compression is 37% 
larger along the c axis than in the basal plane the pressure- induced increase in 
Tc for optimally doped material is thought to occur because pressure moves a band 
associated with the HgOx layer across the Fermi energy, creating new carriers and 
"metallizing" the blocking layer I^S], Even when such dramatic effects do not 
occur, the anisotropic compression in layered materials, such as the copper oxides, 
can move critical features of the density of states (such as the van Hove singularity) 
with respect to the Fermi energy resulting in changes in the carrier density ||30[] . 

MgB2 presents a situation where similar phenomena could occur. The supercon- 
ductivity is thought to result from strong electron-phonon coupling to a particular 
feature of the electronic structure associated with boron a bonds which lies close 



to the Fermi energy [0. The layered structure of MgB2, characterized by Mg-B 
bonds along the c axis and B-B bonds in the basal plane, is expected to compress 
anisotropically. Thus, accurate structural data versus pressure are needed to evalu- 
ate the pressure-induced changes in the electronic structure, as well as the changes 
in phonon frequencies and electron-phonon coupling, and how these might contribute 
to the pressure dependence of Tc. 

Within months after the discovery of superconductivity in MgB2, several groups 
reported structural measurements versus pressure [^, ^ 0. The compression is 
clearly anisotropic, but quantitative agreement among the experimental results for the 
bulk modulus and compression anisotropy was poor. Measurements made in helium 
gas appear to exhibit the largest compression anistropy {[{dc/dP)/co]/{da/dP)/ao]} 
= 1.64(4) and 1.9(3) while measurements made in other fluids yield lower 



values ~ 1.5 in a methanol:ethanol:water mixture |31| and ~ 1.4 in silicone oil ||32|| . 
Some of these differences may be due to the degree to which the pressure fluid is truely 
hydrostatic. Errors in the accurate in situ measurement of lattice parameters and the 
extrapolation of the results to zero pressure could also contribute to the differences. 
Some authors |2^, ^ have concluded that different samples of MgB2 can exhibit 
different pressure- dependent behavior. It has been speculated that samples may 
differ in the amount of Mg or B vacancies, but there is no clear evidence that such 
deviations in stoichiometry are possible in MgB2. An alternative is that impurity 
phases such as MgB4, or elemental Mg or B, distributed at grain boundaries or at 
the center of grains, modify the pressure, and amount of shear, seen by individual 
crystallites of MgB2 in sintered grains when pressure is applied. 

In this paper, we report parallel in situ neutron powder diffraction and Tc{P) 
measurements versus pressure on the same MgB2 sample in a He-gas apparatus to 
0.6 GPa, thus avoiding any problems with non-hydrostatic pressure fluids or sample 
dependent differences. In addition, we present measurements in a helium-loaded 
diamond-anvil-cell to 20 GPa on the same sample. The high precision achieved in 
these measurements allows a quantitative interpretation of the change in Tc versus 
the changes in structure. 
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2 EXPERIMENTAL METHODS 



2.1 Sample Preparation 

The powder sample of MgB2 for these studies was made using isotopically-enriched 
^^B (Eagle Picher, 98.46 atomic % enrichment). A mixture of ^^B powder (less than 
200 mesh particle size) and chunks of Mg metal was reacted for 1.5 hours in a capped 
BN crucible at 800°C under an argon atmosphere of 50 bar. As discussed below, the 
resulting sample displays sharp superconducting transitions in the ac susceptibility 
with full shielding. At ambient pressure the temperatures of the superconducting 
onset and midpoint lie at 39.25 K and 39.10 K, respectively. Since this sample 
contains isotopically pure ^^B, a temperature shift of ATc ~ 0.2 K should be added 
to our Tc values before comparing them with those from other groups using samples 
not isotopically enriched (^°-^^B). 



2.2 Neutron Powder Diffraction Measurements 



Neutron powder diffraction measurements were made on the Special Environment 
Powder Diffractometer at the Intense Pulsed Neutron Source, Argonne National Lab- 
oratory in a He-gas pressure cell |3^ at room temperature. Typical data col- 
lection times were one hour at each pressure. Pressures were measured continuously 
at the pumping station, connected to the pressure cell by a capillary line, and are 
accurate and stable within 0.02 GPa. The data were analyzed by the Rietveld tech- 
nique using the GSAS code In initial refinements, the Mg/^^B ratio was refined. 
There was no indication of non stoichiometry within a refinement precision of about 
0.5%. Fig. 1 shows the raw data and refined diffraction pattern at 0.63(2) GPa. 
The sample is single phase and the diffraction pattern is nicely fit with peak widths 
near the instrumental resolution. This is true at all pressures. There is no evidence 
for any structural transitions. Tc(P). 



2.3 Measurements in He-Gas Apparatus 

The measurements of Tc{P) to 0.7 GPa were carried out using a He-gas high-pressure 
system (Harwood). The pressure is determined by a calibrated manganin gauge at 
room temperature (RT) located in the compressor system. The CuBe pressure cell 
(Unipress) is inserted into a closed-cycle cryocooler (Leybold) with a base temperature 
of 2 K and connected to the compressor system by a 3 mm O.D. x 0.3 mm I.D. CuBe 
capillary tube approximately 3 m long. To minimize shear stresses on the sample 
when the helium pressure medium freezes, a technique developed by Schirber |3^ 



is applied, whereby the top of the 15 cm long pressure cell and the capillary tube 
are kept at a slightly higher temperature than the bottom so that helium freezes 
from the bottom up around the sample as the pressure cell is slowly cooled (30 min) 
through the melting curve of helium. The pressure in the cell can be changed at 
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any temperature above the melting curve Tm{P) of the hehum pressure medium (for 
example, ^ 13.6 K at 0.1 GPa and ^ 38.6 K at 0.50 GPa [|0|). For pressures 
P > 0.5 GPa, Tm lies above the superconducting transition temperature of MgB2 and 
the sample is in frozen helium during the Tc measurement; the slight pressure drop 
(few 0.01 GPa's) on cooling in the solid helium pressure medium from Tm to Tc is 
estimated using the known isochores of He All pressures are determined at T^. 

The superconducting transition of the 8.12 mg MgB2 powder sample is measured 
by the ac susceptibility technique using a miniature primary /secondary coil system 
located inside the 7 mm I.D. bore of the pressure cell. An EG&G 5210 lock-in 
amplifier with a transformer preamplifier is used at 0.113 Oe (rms) field and 1,023 
Hz. A small Pb sphere with 1.76 mm dia (38.58 mg) is also inserted in the coil 
system for susceptibility calibration purposes; for selected data the superconducting 
transition temperature of this Pb sphere is used as an internal manometer |^| to 
check the pressure indicated by the external manganin gauge. 



Fig. 1. Observed time-of-flight neutron powder diffraction data and best-fit Ri- 
etveld refinement profile for MgB2 at 0.63(2) GPa. Data collection time was one 
hour. Crosses (+) are the raw data. The solid line is the calculated profile. Tick 
marks indicate the positions of all allowed reflections. A difference curve (observed 
minus calculated) is plotted at the bottom. 



2.4 Tc{P) Measurements in Diamond- Anvil- Cell 

Tc{P) can be determined to much higher pressures using a helium- loaded diamond- 
anvil-cell made of hardened Cu-Be alloy fitted with 1/6-carat diamond anvils and 0.5 
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mm culet diameter. The MgB2 sample (80 x 80 x 25 yum^) together with several small 
ruby spheres (5-10 fim dia.) are placed in a 240 /xm dia. hole drilled through 
the center of the TaW gasket. The pressure in the gasket hole can be changed at 
any temperature from 1.6 K to RT. Temperature is measured by calibrated Pt and 
Ge thermometers thermally anchored to the top diamond. The pressure in the cell 
can be determined at any temperature below room temperature (RT) to within 0.2 
GPa by measuring the pressure-induced shift in the ruby Rl fluorescence line. The 
pressure is normally measured at temperatures close to the Tc of MgB2. 

The superconducting transition itself is determined inductively to ± 0.1 K us- 
ing two balanced primary/ secondary coil systems connected to a Stanford Research 
SR830 digital lock-in amplifier. The ac susceptibility studies were carried out using 
a 3 G (r.m.s.) magnetic field at 1000 Hz. Over the transition the signal changed by 
~ 3 nV with a background noise level of ~ 0.2 nV. Further details of the He-gas and 



diamond-anvil-cell high-pressure techniques are given elsewhere ^3], P 



3 RESULTS OF EXPERIMENT 

3.1 Pressure-Dependent Structural Properties 

The simple hexagonal structure of MgB2 (space group P6/mmm, No. 191) is shown 
in Fig. 2. The structure contains graphite-like boron layers which are separated by 
hexagonal close-packed layers of metals. The center of a hexagonal boron ring lies 
both directly above and below each metal. 

The variation of the a and c lattice parameters vs. pressure is shown in Fig. 3. 
Over the pressure range of this study, the changes are linear and can be expressed as 

a = ao[l - 0.00187(4)P] and c = Co[l - 0.00307(4)P], (1) 

where Oq = 3.08489(3) and Cq = 3.52107(5) are the zero-pressure lattice parameters 
and P is the pressure in GPa. Numbers in parenthesis are standard deviations 
of the last significant digit. The bulk modulus [VoldP / dV)] obtained from these 
measurements is 147.2(7) GPa. 

Loa and Syassen used electronic structure calculations vs. cell volume to calculate 
a bulk modulus of 140.1(6), in good agreement with the experimental result. They 
also calculated the pressure dependence of the c/a ratio, getting a result in nice 
agreement with the observed compression anisotropy. 

The compression anisotropy, defined as [{dc/dP)/co]/ {da/dP)/ao], is 1.64(4). Com- 
pression along the c axis is 64% larger than along the a axis, consistent with the 
comparatively weaker (Mg-B) bonds that determine the c axis length. A similar 
anisotropy, but not as large, has been reported in the refractory diboride TiB2 PB| , 



which is of considerable technological interest because of its high elastic moduli, high 
hardness, and high electric conductivity. By comparison, the compression anisotropy 
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in the layered cuprate YBa2Cu307 is about a factor of two Not surprisingly, 

the intrinsic compression anisotropy is not observed when pressure measurements are 
made in non-hydrostatic media. Recent room-temperature x-ray diffraction measure- 



ments in diamond anvil cells using methanol:ethanol:water |^T|] and silicone oil as 
the pressure fluids gave anisotropies of 1.5 and 1.4, respectively. An x-ray diffraction 
study to much higher pressures using helium as the pressure fluid in a diamond anvil 
cell gives a compression anisotropy of 1.9(3), in agreement with our result within 
the error bars. 



Fig. 2. Crystal structure of MgB2 [AlB2-type structure; hexagonal space group 
P6/mmm, No. 191, with Mg at (0, 0, 0) and B at (1/3, 2/3, 1/2)] viewed along the 
c axis (top) and perpendicular to an a axis (bottom). Small spheres are B atoms; 
larger spheres are Mg atoms. 
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Fig. 3. Normalized a and c lattice parameters vs. pressure at room temperature 



for MgB2 based on neutron diffraction measurements from Ref. |Q at five pressures 
using helium as the pressure transmitting medium. Standard deviations of the 
individual points are smaller than the symbols. The straight lines are linear least- 
sqaures fits to the data. 

3.2 Pressure-Dependent Superconducting Properties 
3.2.1 Tf.{P) Measurements in the He-Gas System 

In Fig. 4 we show representative examples of the superconducting transition for 
MgB2 in the ac susceptibihty at both ambient and high pressure in the He-gas system 
With increasing pressure the narrow transition is seen to shift bodily to lower 
temperatures, allowing a determination of the pressure-induced shift in to within 
± 10 mK. Remarkably, close inspection of the data for 0.50 GPa reveals a slight 
jog in the transition curve near its midpoint, accurately marking the position of the 
melting curve of helium (T^ — 38.6 K) at this pressure. 
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Fig. 4. Real part of the ac susceptibility of MgB2 versus temperature at ambient 
and high pressures from Ref. |^^. The applied magnetic field is 0.113 Oe (rms) 
at 1,023 Hz. Intercept of straight tangent lines defines superconducting onset at 
ambient pressure T°"*^*(0) ~ 39.25 K, with the superconducting midpoint T™'^(0) ~ 
39.10 K. No correction is made for demagnetization effects. 



In Fig. 5, the dependence of on pressure is seen to be highly linear dT^/dP ~ 
—1.11(2) K/GPa. Data were obtained following pressure changes at both RT (un- 
primed data) and low temperature (primed data). The dependence of on pressure 
thus does not depend on the pressure/temperature history of the sample. Such his- 
tory effects are rare in superconductors without pressure-induced phase transitions, 
but do occur in certain high- Tc oxides containing defects with appreciable mobility 



at RT 48 



In selected loosely bound solids with large molecular units, like Cgo, helium atoms 
are able to intercalate inside when pressure is applied, diminishing the pressure- 
induced changes in the sample properties [^. An analysis of the MgB2 structure 
readily reveals that its hexagonal unit cell is tightly packed with insufficient space for 
helium atoms to readily travel through. To verify that helium does not intercalate 
inside MgB2 under pressure, we carried out a parallel experiment to 0.077 GPa with 
neon gas instead of helium. In analogy with the results on Cgo the intercalation 
of the larger neon atoms into MgB2 would be more difficult than for helium. The fact 
that the pressure derivative dT^/dP is the same for both helium and neon confirms 
the absence of intercalation effects in the present experiments. 
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Fig. 5. Superconducting transition temperature onset versus applied pressure from 
Ref. Numbers give order of measurement. Data for pts. 2', 6, 8', and 11 are 



shown in Fig. 4. A typical error bar for (±0.01 K) is given in lower left corner; 
the error in pressure is less than the symbol size. Pressure was either changed at RT 
(unprimed numbers) or at low temperatures ~ 60 K (primed numbers). 

Since for pressures less than 0.5 GPa the sample is surrounded by fluid helium 
during the measurement, the measured slope dT^/dP ~ -1.11 K/GPa to this 
pressure gives the true hydrostatic pressure dependence for MgB2. For P > 0.5 GPa 
the sample is in frozen helium at temperatures near Tc, but, as seen in Fig. 5, no 
change in the pressure dependence Tc(P) is observed. This is not surprising since 
solid helium is the softest solid known; in addition, the shear stresses are held to a 



minimum by the carefully controlled manner |^ in which solid helium is allowed to 
freeze around the sample. 

A similar pressure derivative dT^/dP ~ —1.07 K/GPa to ours has very recently 
been obtained by Lorenz et al. in He-gas studies to 0.8 GPa on a MgB2 sample 



synthesized to stoichiometry with superconducting midpoint at T™'^(0) ~ 39.2 K. 
These authors also reexamined in a He-gas system the same sample studied earlier 
22[ with T™'^(0) ~ 37.5 K and find dTjdP ~ -1.45 K/GPa which they report 
agrees within experimental error with their previous result dTc/dP ~ —1.6 K/GPa 
to 1.8 GPa in a piston-cylinder device with Fluorinert FC77 pressure medium. This 
appears to imply that the shear stresses from frozen Fluorinert have little effect on 
the pressure dependence of Tc in the pressure range to 1.8 GPa. 
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Choi et al. have recently carried out resistivity studies to 1.5 GPa pressure in 
daphne-kerosene pressure medium, obtaining dT^/dP ^ —1.36 K/GPa. The results 
of all known Tc{P) measurements on MgB2 are summarized in the Table. 



3.2.2 Tc{P) Measurements in the Diamond-Anvil System 

In Fig. 6 we show the dependence of on pressure to 20 GPa for MgB2 using 
a diamond-anvil-cell with dense helium pressure medium [Q, thus extending the 
pressure range of the above He-gas studies nearly thirtyfold. Tc is seen to decrease 
nearly linearly with pressure to 10 GPa, consistent with the rate -1.11 K/GPa (dashed 
line), but begins to display a positive (upward) curvature at higher pressures. As 
will be discussed below, this deviation originates from the increasing lattice stiffness 
of MgB2 at higher pressure. In these experiments the pressure was always changed 
at RT, but measured at temperatures near Tc. 



Fig. 6. Superconducting transition temperature midpoint T™'' versus pressure to 20 
GPa from diamond-anvil-cell measurements in Ref. Data with filled circles (•) 

taken for monotonically increasing pressure, with open circles (o) for monotonically 
decreasing pressure. The straight dashed line has slope -1.11 K/GPa. 



In Fig. 6 it is seen that the width of the superconducting transition gradually 
increases from ~ 0.3 K for P < 10 GPa to 0.9 K at 19.2 GPa, increasing somewhat 
further for the data with decreasing pressure. This increase in width ATc is seen to 
be usually accompanied by a slight broadening of the ruby Rl fluorescence line; both 
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broadening effects point to a pressure gradient of approximately ± 0.3 GPa (± 1.5%) 
at the highest pressures. The magnitude of the shear stresses on the sample would 
be expected to be larger in the diamond-anvil-cell than in the He-gas experiment 
since the pressure range is much greater; above 12 GPa helium freezes at RT so the 
diamonds must push on solid helium to increase the pressure further. In addition, in 
the diamond-anvil-cell it is not possible to cool slowly through the melting curve of 
helium with a well-defined temperature gradient. However, the data in Fig. 6 give 
no clear indication for shear stress effects on Tc at any pressure. 

Very recently Tissen et al. |]^ have carried out ac susceptibility measurements 
in a diamond-anvil-cell to 28 GPa on a MgB2 sample with T^'"*'^(0) ~ 37.3 K at 
ambient pressure. They find an initial slope dT^/dP ~ —2 K/GPa, decreasing 
to 11 K at 20 GPa and 6 K at 28 GPa, a 50% greater decrease than observed by 
either us (see Fig. 6) or Monteverde et al. [24|. They also report that the pressure 
dependence Tc{P) shows a bump near 9 GPa which they speculate may arise from 
an electronic Lifshitz transition. We suggest that shear stress effects may also play a 
role in their measurements. At 20 GPa the width in their superconducting transition 
has increased by ~ 3 K which would correspond to a pressure gradient of ~ 3.5 
GPa, an order of magnitude higher than in our helium-loaded diamond-anvil-cell 
measurements. 

The degree to which shear stresses affect the data of Monteverde et al. is 
unknown. However, since shear stresses are potentially much larger in solid pressure 
media such as steatite than in frozen fluids such as helium or Fluorinert, it would 
seem likely that they are responsible for at least part of the widely differing Tc{P) 
dependences to 25 GPa observed in three of their four experiments. As discussed in 
the Introduction, the initial pressure derivative dTc/dP in the present experiment 
(-1.11 K/GPa) differs significantly from those (-0.3 to -2.0 K/GPa) obtained by other 
groups using pressure media which are either solid at RT or readily freeze upon 



cooling [El, E2|, E3|, pO, E^. It is not yet clear whether these widely varying results 



reflect differences in the make-up of the samples or differing degrees of shear stress 
exerted on the samples by the various frozen or solid pressure media. An inspection 



of the data in the Table suggests a possible correlation voiced by Tissen et al. |5T 
that larger values of \dTc/dP\ are associated with lower ambient-pressure values of T^. 
However, it is difficult to accurately compare values determined in ac susceptibility 
and electrical resistivity measurements, the latter usually lying higher; in addition, in 
the ac susceptibility the value of Tc may depend somewhat on the ac field strength. 
Further experimentation under carefully controlled conditions is clearly necessary to 
investigate this possible correlation. 
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Table. Summary of available high-pressure Tc{P) data on MgB2. values are at 
ambient pressure from superconducting midpoint in ac susceptibility Xac and electrical 
resistivity p measurements. dTc/dP is initial pressure derivative. P™'^^(GPa) is the 
maximum pressure reached in experiment . 
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4 DISCUSSION 

The present studies of both the superconducting and structural properties of MgB2 
under hydrostatic pressure were carried out on the same high quality MgB2 sample 
used in the He-gas measurements to 0.7 GPa. These combined studies thus allow an 
accurate determination of the change in Tc with unit cell volume V for comparison 
with theory. The change in Tc with V is given by 



dhiTr B /dT, 



dlnV T, V dP 



+4.16(8), (2) 



using the above values dTjdP ~ -1.11(2) K/GPa, B = 147.2(7) GPa, and = 
39.25 K. This value of dlnTc/dlnV is somewhat smaller than that (+6.6) obtained 
by Neaton and Perali in an estimate based on density functional theory. 



We will now discuss the implications of this result for the nature of the supercon- 



ducting state in MgB2. First consider the McMillan equation ^4 



1.20 ^\A-/i*(l + 0.62A) j ' ^ ^ 

valid for strong coupling (A < 1.5), which connects the value of Tc with the electron- 
phonon coupling parameter A, an average phonon frequency (uj) , and the Coulomb 
repulsion /x*, which we assume to be pressure independent [Q. The coupling param- 
eter is defined by A = N{Ef) {P) /[M(co'^)], where N{Ef) is the electronic density 
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of states at the Fermi energy, {P) the average squared electronic matrix element, 
M the molecular mass, and {uo"^) the average squared phonon frequency. Taking the 
logarithmic volume derivative of Tc in Eq. (3), we obtain the simple relation 

-7 + A <^ TTTT + 27 L (4) 



d\nV ' \dhiV 

where 7 = —din (uj) /d\nV is the Griineisen parameter, rj = N{Ef) (P) is the Hop- 
field parameter |5^, and 

^ ^ 1.04A[l + 0.38/i1 
[A-^*(1 + 0.62A)]'' 

Eq. (4) has a simple interpretation. The first term on the right, which comes from the 
prefactor to the exponent in the above McMillan expression for Tc, is usually small 
relative to the second term, as will be shown below. The sign of the logarithmic 
derivative c? In T^/c? In V^, therefore, is determined by the relative magnitude of the 
two terms in the curly brackets. 

The second "electronic" term in Eq. (4) involves the logarithmic volume derivative 
of the Hopfield parameter t] = N{Ef) (/^), an "atomic" property which can be calcu- 
lated directly in band-structure theory |Q. In his landmark paper [^, McMillan 
demonstrated that whereas N{Ef) and (/^) individually may fluctuate appreciably 
as one element is substituted for another across a transition-metal alloy series and 
the d-electron count varies, their product rj = N{Ef) {P) changes only gradually, i.e. 
rj is a. well behaved "atomic" property. One would thus anticipate that rj changes 
in a relatively well defined manner under pressure, reflecting the character of the 
electrons near the Fermi energy. An examination of the body of high-pressure data 
on simple s,p-metal superconductors, in fact, reveals that rj normally increases under 
pressure at a rate given hy dlnrj/ d\nV ~ —1 |^. For transition-metal (d-electron) 
superconductors, on the other hand, Hopfield has pointed out that dXnrj/ d\n.V ~ —3 
to -4 pf. 

The second "lattice" term in the curly brackets in Eq. (4) is positive, typically 
27 ^ 3 — 5. Since in simple metal superconductors, like Al, In, Sn, and Pb, this 
positive "lattice" term dominates over the electronic term dhirj/ d\aV ~ —1, and A 
is always positive, the sign of d In Td dlnV is the same as that in the curly brackets, 
namely positive; this accounts for the universal decrease of Tc with pressure due to 
lattice stiffening in simple metals. In selected transition metals the electronic term 
may become larger than the lattice term, in which case d\n.Tc/d\n.V is negative 
and Tc would be expected to increase with pressure, as observed, for example, in 
experiments on V and La 

Let us now apply Eq. (4) in more detail to a canonical BCS simple-metal super- 
conductor. In Sn, for example, Tc decreases under pressure at the rate dTc/dP ~ 



-0.482 K/GPa which leads to dlnTjdlnV ~ +7.2 ETI. We note that this value of 
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dlnTc/dlnV is almost twice as large as that for MgB2 (see Eq. (2)); this is exactly 
what is expected from Eq. (4) since A increases for decreasing values of T^. Inserting 
for Sn Tc(0) ~ 3.73 K, (a;) ~ 110 K Q, and /i* = 0.1 into the above McMillan 
equation, we obtain A ~ 0.69 from which follows that A ~ 2.47. Inserting the above 
values into Eq. (4) and setting dlnrj/dlnV ~ —1 for simple metals, we can solve Eq. 
(4) for the Griineisen parameter to obtain 7 ^ +2.46, in reasonable agreement with 
experiment for Sn (7 ^ +2.1) |^T[]. Similar results are obtained for other conventional 
simple metal BCS superconductors. 

We now repeat the same calculation with the McMillan equation for MgB2 using 



the logarithmically averaged phonon energy from inelastic neutron studies |TT| {cu) 



670 K, Tc(0) ~ 39.25 K, and fx* = 0.1, yielding A ~ 0.90 and A ~ 1.75 from Eqs. 
(3) and (5), respectively. Our estimate of A ^ 0.90 agrees well with those of other 
authors [?, |16|. Since the pairing electrons in MgB2 are believed to be s,p in character 



15] , p!7| , ?, Q , we set dlnrj / d\nV —1, a value close to dlnrj/dlnV = Bd In 17/ dP ~ 
-0.81, where B = 147.2 GPa from Ref. fSSf and d\nr]/dP ^ +0.55 %/GPa from 
first-principles electronic structure calculations by Medvedera et al. |^T|]. Inserting 
the values of dhiTc/d\n.V = +4.16, A = 1.75, and d\nr)/dlnV = —1 into Eq. 
(4), we find 7 ~ 2.36, in reasonable agreement with the value 7 ~ 2.9 from Raman 
spectroscopy studies |Q or 7 2.3 from ab initio electronic structure calculations 
on MgBa 



In spite of the significant compression anisotropy, electronic structure calculations 
based on the high-pressure structural data show that the electronic structure does 
not change much at high pressure the calculations show that the electric field 
gradient in MgB2 is essentially independent of pressure up to 10 GPa. As the electric 
field gradient is a very sensitive characteristic of the electronic charge distribution, one 
may conclude that no large changes in the partial charges of the B 2p states and boron 
electronic structure take place under pressure. Further results from theory support 
this conclusion. Medvedera et al. [^] find the Hopfield parameter for MgB2 to only 
depend weakly on pressure dlnrj/dP ~ +0.55 %/GPa. The change in the electronic 
density of states d\nN{Ef)/dP is also estimated to be very small: Loa and Syassen 
PI (-0.31 %/GPa), Medvedera et al. §^ (-0.51 %/GPa), and Vogt et al. (-0.38 
%/GPa). Assuming B = 147.2 GPa, one thus obtains d\nN{Ef)/d\nV ~ +0.46, 
+0.75, and +0.56, respectively. These values are near that (+0.67) expected for a 
3D free electron gas. Since dlnrj/dP = dlnN{Ef)/dP + din (/^) /dP, these results 
imply that the average squared electronic matrix element (J^) in MgB2 increases under 
pressure at the approximate rate of only +1 %/GPa. The sign and magnitude of the 
changes in N{Ef) and (/^) under hydrostatic pressure for MgB2 are comparable to 
those found for simple s,p-metal superconductors. Larger changes are anticipated if 
uniaxial pressure is applied . The main reason for the observed decrease of Tc with 
pressure is not an electronic effect, but a strong pressure enhancement of the phonon 
frequencies, an effect which has been directly observed in Raman measurements p4 



Taken as a whole, the above results thus give considerable evidence that the super- 
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conducting state of MgB2 is strongly related to that in simple s,p-metal superconduc- 
tors like Al, Sn, In, and Pb which exhibit BCS phonon-mediated superconductivity. 
This is not to say that superconductivity in MgB2 is identical to that in the simple 
metals. Extensive specific heat [|10| and high-resolution photoemission studies 
on MgB2 give evidence for a multicomponent superconducting gap. 

The above analysis is based on the results of the present high-pressure studies 
using the He-gas technique to 0.7 GPa . We now consider the diamond-anvil-cell 
data to 20 GPa in Fig. 6. For comparison to theory it is advantageous to use the 
Murnaghan equation-of-state to convert pressure to relative volume V jV^ 



B'P 



-1 -l/B' 



(6) 



where we use the value B = 147.2 GPa from Ref. |Q and the canonical value 
B' = dB/dP = 4 supported by a recent calculation [^. In Fig. 7 we replot the 
data from Fig. 6 as Tc versus relative volume V/Vq. The maximum pressure applied 
in the present experiment (19.2 GPa) results in a volume decrease of ~ 10%. Much 
of the nonlinearity in the versus pressure plot in Fig. 6 appears to disappear when 
Tc is plotted versus V/Vq. 

We now compare the versus V/Vq dependence in Fig. 7 to the result from the 
He-gas data which yields the initial volume dependence dlnTc/dlnV ~ +4.16 given 
in Eq. (2). If we assume this relation holds at all pressures, then we can integrate it 
to obtain 



(39.25 K) V^c 



(7) 



which is plotted as the upper solid line in Fig. 7. This volume dependence must 
be accurate for small pressures where V/Vq ~ 1, corresponding to the pressure de- 
pendence dTc/dP = —1.11 K/GPa from the He-gas data, but rises well above the 
experimental data at higher pressures. 

Another way to extrapolate the He-gas data to higher pressures is to assume that 
Tc varies linearly with volume change AV, yielding from Eq. (7) 

' ° ^ - ■ 4.16-— = -3.16 + 4.16— , (8) 



39.25 K V K) / V K) / V 

which is plotted as the straight dashed line in Fig. 7. As they must, the upper solid 
and dashed lines agree exactly near V jVi^ = 1. The dashed fine is seen to lie above 
the experimental data points at higher pressures and to extrapolate to Tc = K for 
V/Vo = 0.76 which corresponds to an applied pressure of ~ 75 GPa. A least-squares 
straight line fit through all data in Fig. 7 leads to the estimate that Tc = K for 
P ^ 60 GPa. 
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Fig. 7. Tc data from figure 6 plotted versus relative volume V/Vq (from Ref. [Q). 
See text for explanation of solid and dashed lines. 



It is not surprising that is not a linear function of V/Vq to very high pressures. 
In the McMillan formula in Eq. (3) Tc depends exponentially on the solid state pa- 
rameters and it is the relatively small changes in these parameters which lead to the 
large change in under pressure see in Figs. 6 and 7. As pointed out by Chen 
et al. ||55|, a more appropriate method to estimate the dependence of Tc on relative 
volume V/Vq would thus be to integrate the volume derivatives of these parame- 
ters 7 = -dln(cj) /dhiV = +2.36, d\nX/d\nV = d\nT]/d\nV - d\n{uj^) /d\nV = 
-1 - 2(-2.36) = +3.72 to obtain {lu) = (670 K){V/Vo)-^-^^ and A = 0.90(l^/K))^-^^. 
Inserting these two volume dependences in the McMillan equation, and assuming 
fi* = 0.1 is independent of pressure [^, we obtain the dependence of on relative 
volume shown as the lower solid line in Fig. 7. The agreement with the experimental 
data is quite impressive. Note that according to this estimate approximately 50 GPa 
pressure would be required to drive Tc to below 4 K. A similar calculation was very 
recently carried out by Chen et al. over a much wider pressure range; this paper 
also contains a detailed discussion of the pressures dependences of i], A, and /i*. The 
good agreement between the experimental data to 20 GPa and the predictions of the 
McMillan formula using the volume dependences determined from the He-gas high- 
pressure data to 0.7 GPa lends additional evidence that superconductivity in MgB2 
originates from standard BCS phonon-mediated electron pairing. 

In conventional metals, electron-phonon scattering makes the dominant contri- 
bution to the temperature-dependent electrical resistivity p{T). At sufficiently high 
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temperatures, Bloch-Griineisen theory gives a linear dependence on temperature 
Prt = bT, where h oc r~^6^^, is the radius of the Wigner-Seitz sphere, and 6d 
is the Debye temperature. Near RT Choi et al. find the electrical resistivity of 



MgB2 to increase linearly with temperature. Under pressure these authors find that 
the RT electrical resistivity decreases under pressure at the rate dhi prt / dP ~ —3 
%/GPa. Using the bulk modulus B = 147.2 GPa, this yields dlnpRx/dlnV ^ 
+4.42. Taking the logarithmic volume derivative of the above Bloch-Griineisen ex- 
pression and using the free-electron expression for and setting 7 = 2.36 from above, 
we obtain din pjiT/dlnV = 27 — 2/3 = +4.05, in surprisingly good agreement with 
the measured value. It is significant that the same value of the Griineisen parameter 
yields the pressure dependence of the electron-phonon interaction which accounts for 
both T,(P) and pRriP)- 

At first glance the present results appear to be inconsistent with the hole super- 
conductivity model of Hirsch and Marsiglio [0, ^ which predicts that Tc should 
increase with pressure if there is no change in the doping level of holes. Indeed, the 
pressure-induced change in the concentration of hole-carriers in the boron cr-band is 



estimated to be extemely small |^, Further experiments, such as high-pressure 
Hall effect measurements, are necessary to determine what, if any, change in the 
carrier concentration occurs. The success of the above analysis of the dependence 
of Tc on pressure gives further evidence that MgB2 is an extraordinary supercon- 
ductor which makes the most out of its conventional BCS electron-phonon pairing 
interaction. 



Acknowledgments 

The authors are grateful to S. Short for assistance with the preparation of figures 
for the paper and to X. J. Chen for providing a preprint of his recent paper. The 
authors would like to thank S. Klotz for providing the ruby spheres. Work at Wash- 
ington University supported by NSF grant DMR-0101809 and that at the Argonne 
National Laboratory by the U.S. Department of Energy, Office of Science, contract 
No. W-31-109-ENG-38. 



REFERENCES 



References 

[1] J. G. Bednorz and K. A. Mliller, Z. Physik, 64 (1986) 189. 

[2] A. F. Hebard, M. J. Rosseinsky, R. C. Haddon, D. W. Murphys, S. H. Glarum, 
T. T. M. Palstra, A. P. Ramirez, and A. R. Kortan, Nature, 350 (1991). 

[3] J. H. Sch5n, Ch. Kloc, and B. Batlogg, Nature, 406 (2000) 704; ibid., Nature, 408 
(2000) 549. 



18 



[4] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and J. Akimitsu, Nature, 
410 (2001) 63. 

[5] D. C. Larbalestier, L. D. Cooley, M. O. Rikel, A. A. Polyanskii, J. Jiang, S. 
Patnaik, X. Y. Cai, D. M. Feldmann, A. gurevich, A. A. Squitieri, M. T. Naus, C. 
B. Eom, E. E. Hellstrom, R. J. Cava, K. A. Regan, N. Rogado, M. A. Hayward, 
T. He, J. S. Slusky, R Khalifah, K. Inumaru, M. Haas, Nature, 410 (2001) 186. 

[6] P. C. Canfield, D. K. Finnemore, S. L. Bud'ko, J. E. Ostenson, G. Lapertot, C. 
E. Cunningham, and C. Petrovic, Phys. Rev. Lett., 86 (2001) 2423. 

[7] S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunningham, N. Anderson, P. C. 
Canfield, Phys. Rev. Lett., 86 (2001) 1877. 

[8] D. G. Hinks, H. Glaus, J. Jorgensen, Nature, 411 (2001) 457. 

[9] R. K. Kremer, B. J. Gibson, K. Ahn, preprint |cond-mat / 0102432] . 

[10] Y. Wang, T. Plackowski, and A. Junod, Physica C, 355 (2001) 179. 

[11] R. Osborn, E. A. Goremychkin, A. I. Kolesnikov, D. G. Hinks, Phys. Rev. Lett., 
87 (2001) 017005. 

[12] T. Yildirim, O. Giilseren, J. W. Lynn, C. M. Brown, T. J. Udovic, H. Z. Qing, 
N. Rogado, K. A. Regan, M. A. Hayward, J. S. Slusky, T. He, M. K. Haas, P. 
Khalifah, K. Inumaru, R. J. Cava, Phys. Rev. Lett., 87 (2001) 037001. 

[13] H. Kotegawa, K. Ishida, Y. Kitaoka, T. Muranaka, and J. Akimitsu, Phys. Rev. 
Lett., 87 (2001) 127001. 

[14] T. Takahashi, T. Sato, S. Souma, T. Muranaka, J. Akimitsu, Phys. Rev. Lett., 
86 (2001) 4915. 

[15] J. Kortus, I. I. Mazin, K. D. Belashchenko, V. P. Antropov, L. L. Boyer, Phys. 
Rev. Lett., 86 (2001) 4656. 

[16] J. M. An and W. E. Pickett, Phys. Rev. Lett., 86 (2001) 4366. 

[17] N. I. Medvedeva, A. L. Ivanovskii, J. E. Medvedeva, A. J. Freeman, Phys. Rev. 
B, 64 (2001) 020501. 

[18] O. F. de Lima, R. A. Ribeiro, M. A. Avila, C. A. Cardoso, and A. A. Coelho, 
Phys. Rev. Lett., 86 (2001) 5974; C. U. Jung, J.-H. Choi, P. Chowdhury, K. H. 
P. Kim, M.-S. Park, H.-J. Kim, J. Y. Kim, Z. Du, M.-S. Kim, W. N. Kang, S.-Ik 
Lee, G. Y. Sund, and J. Y. Lee, 

[19] S. Lee, H. Mori, T. Masui, Yu. Eltsev, A. Yamamoto, and S. Tajima, preprint 
|cond-mat/0 105545 . 



19 



[20] See, for example: U. Welp, W. K. Kwok, G. W. Crabtree, K. G. Vandervoort, 
and J. Z. Liv, Phys. Rev. Lett., 62 (1989) 1908. 

[21] See, for example, P. E. Seiden, Phys. Rev., 179 (1969) 458. 

[22] B. Lorenz, R. L. Meng, C. W. Chu, Phys. Rev. B, 64 (2001) 012507. 

[23] E. Saito, T. Taknenobu, T. Ito, Y. Iwasa, K. Prassides, and T. Arima, J. Phys.: 
Condens. Matter, 13 (2001) L267. 

[24] M. Monteverde, M. Niifiez-Regueiro, N. Rogado, K. A. Regan, M. A. Hayward, 
T. He, S. M. Loureiro, R. J. Cava, Science, 292 (2001) 75. 

[25] S. Klotz and J. S. Schilling, Physica C, 209 (1993) 499; J. S. Schilling and S. 
Klotz, in: Physical Properties of High Temperature Superconductors, Vol. Ill, ed. 
D. M. Ginsberg (World Scientific, Singapore, 1992) p. 59. 

[26] T. Ishiguro, H. Ito, Y. Yamauchi, E. Ohmichi, M. Kubota, H. Yamochi, G. Saito, 
M. V. Kartsovnik, M. A. Tanatar, Yu. V. Sushko, G. Yu. Logvenov, Synthetic 
Metals, 85 (1997) 1471; D. Jerome and H. J. Schulz, Adv. Phys., 31 (1982) 299. 

[27] B. A. Hunter, J. D. Jorgensen, J. L. Wagner, P. G. Radaelh, D. G. Hinks, H. 
Shaked, R. L. Hitterman, R. B. VonDreele, Physica C, 221 (1994) 1. 

[28] D. L. Novikov, O. N. Myrasov, A. J. Freeman, Physica C 222, 38 (1994); D. J. 
Singh & W. E. Pickett, Physica C, 233 (1994) 237. 

[29] J. D. Jorgensen, D. G. Hinks, O. Chmaissem, D. N. Argyriou, J. F. Mitchell, 
B. Dabrowski, in: Recent Developments in High Temperature Superconductivity, 
edited by J. Klamut, B. W. Veal, B. M. Dabrowski, P. W. Klamut, (Springer- 
Verlag, Berlin, 1996) p. 1. 

[30] D. L. Novikov & A. J. Freeman, in Recent Developments in High Temperature 
Superconductivity, edited by J. Klamut, B. W. Veal, B. M. Dabrowski, P. W. 
Klamut, (Springer- Verlag, Berlin, 1996) p. 17. 

[31] T. Vogt, G. Schneider, J. A. Hriljac, G. Yang, J. S. Abell, Phys. Rev. B, 63 
(2001) 220505(R). 

[32] K. Prassides, Y. Iwasa, T. Ito, Dam H. Chi, K. Uehara, E. Nishibori, M. Takata, 
S. Sakata, Y. Ohishi, O. Shimomura, T. Muranaka, J. Akimitsu, Phys. Rev. B, 
64 (2001) 012509. 

[33] J. D. Jorgensen, D. G. Hinks, S. Short, Phys. Rev. B, 63 (2001) 224522. 

[34] A. F. Goncharov, V. V. Struzhkin, E. Gregoryanz, J. Hu, R. J. Hemley, Ho-k. 
Mao, G. Lapertot, S. L. Bud'ko, P. C. Canfield, preprint |cond-mat/0104042| . 



20 



[35] B. Lorenz, R. L. Meng, and C. W. Chu, preprint |cond-mat / 1 04303 . 

[36] J. D. Jorgensen, J. Faber Jr., J. M. Carpenter, R. K. Crawford, J. R. Haumann, 
R. L. Hitterman, R. Kleb, G. E. Ostrowski, F. J. Rotella, and T. G. Worlton, J. 
Appl. Crystallogr., 22 (1989) 321. 

[37] J. D. Jorgensen, S. Pei, R Lightfoot, D. G. Hinks, B. W. Veal, B. Dabrowski, 

A. R Paulikas, R. Kleb, Physica C, 93 (1990) 93. 

[38] A. C. Larson and R. B. VonDreele (unpublished) 
[39] J. E. Schirber, Cryogenics, 10 (1970) 418. 
[40] I. L. Spain and S. Segall, Cryogenics, 11 (1971) 26. 
[41] A. Filing and J. S. Schilling, J. Phys., F 11 (1981) 623. 

[42] The ruby spheres, which were annealed at 1500°C for one week followed by 
a cooldown to RT over a three-week period, were prepared by J. C. Chervin, 
Physique des Milieux Condenses, Paris, and contained 3600 ppm Cr'^"'". 

[43] J. S. Schilling, J. Diederichs, S. Klotz, R. Sieburger, in: Magnetic susceptibility of 
Superconductors and Other Spin Systems, edited by R. A. Hein, T. L. Francavilla, 
D. H. Liebenberg (Plenum Press, New York, 1991), p. 107. 

[44] S. Deemyad, J. S. Schilling, J. D. Jorgensen, and D. G. Hinks, Physica C, 361 
(2001) 227. 

[45] I. Loa and K. Syassen, Sohd State Commun., 118 (2001) 279. 

[46] P. S. Spoor, J. D. Maynard, M. J. Pan, D. J. Green, J. R. Hellimann, T. Tanaka, 
Appl. Phys. Lett., 70 (1997) 1959. 

[47] T. Tomita, J. J. Hamlin, J. S. Schilling, D. G. Hinks, and J. D. Jorgensen, Phys. 
Rev. B, 64 (2001) 092505. 

[48] See, for example, S. Sadewasser, J. S. Schilling, and A. M. Hermann, Phys. Rev. 

B, 62 (2000) 9155. 

[49] J. F. Schirber, G. H. Kwei, J. D. Jorgensen, R. L. Hitterman, and B. Morosin, 
Phys. Rev. B, 51 (1995) 12014. 

[50] F. S. Choi, W. Kang, J. Y. Kim, M.-S. Park, C. U. Jung, H.-J. Kim, and S.-Ik 
Lee, preprint |cond-mat / 10445^ . 

[51] V. G. Tissen, M. V. Nefedova, N. N. Kolesnikov, and M. P. Kulakov, preprint 
|cond-mat/0105475| . 



21 



[52] J. J. Hamlin, C. Looney, T. Tomita, J. S. Schilling, D. G. Hinks, and J. D. 
Jorgensen (unpublished) . 



[53 
[54; 

[55 
[56 
[57; 

[58 
[59 
[60 
[61 

[62 
[63 

[64; 

[65 
[66 

[67; 



J. B. Neaton and A. Perali, preprint |cond-mat/010409^ . 

W. L. McMillan, Phys. Rev., 167 (1968) 331; P. B. Allen and R. C. Dynes, Phys. 
Rev. B, 12 (1975) 905. 

X. J. Chen, H. Zhang, and H.-U. Habermeier, preprint. 

J. J. Hopfield, Physica, 55 (1971) 41. 

J. S. Schilling and S. Klotz, in: Physical Properties of High Temperature Su- 
perconductors, Vol. Ill, ed. D.M. Ginsberg (World Scientific, Singapore, 1992) p. 
59. 

T. F. Smith and H. L. Luo, J. Phys. Chem. Sohds, 28 (1967) 569. 

T. F. Smith, J. Phys. F, 2 (1972) 946. 

R. C. Dynes, Sohd State Commun., 10 (1972) 615. 

N. I. Medvedeva, A. L. Ivanovskii, J. E. Medvedeva, A. J. Freeman, and D. L. 
Novokov (preprint |cond-mat/0 10434^ ) . 

D. Roundy, H. J. Choi, H. Sun, S. G. Louie, M. L. Cohen, post deadline session 
on MgB2, Monday evening, March 12, 2001, APS March Meeting in Seattle, 
Washington. 

K. Kobayashi and K. Yamamoto, preprint |cond-mat /0105231| ; X. Wan, J. Dong, 



H. Weng, and D. Y. Xing, preprint |cond-mat /01042T6 



S. Tsuda, T. Tokoya, T. Kiss, Y. Takano, K. Tongano, H. Kitou, H. Ihara, and 
S. Shin, preprint |cond-mat/0104489| . 

See, for example, J. M. Ziman, Principles of the Theory of Solids (Cambridge 
U.P., Cambridge, England, 1964) p. 191; see also, J. S. Schilling and W. B. 
Holzapfel, Phys. Rev. B, 8 (1973) 1216. 

J. E. Hirsch, Phys. Lett. A, 282 (2001) 392. 

J. E. Hirsch and F. Marsiglio, Phys. Rev. B, 64 (2001) 144523. 



22 



This figure "Figlnew.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



This figure "Narfigl.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



This figure "Fig2new.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



This figure "Narfig2.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



This figure "Tc_Vs_Press.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



This figure "Tc_Vs_VolChange.png" is available in "png" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



This figure "narfigS.PNG" is available in "PNG" format from: 



http://arXiv.org/ps/cond-mat/0110268vl 



